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We extend the study of a J p = 2 + ,I = § TrAiV quasibound state [Phys. Rev. D 78, 014013 

(2008)] by solving nonrelativistic Faddeev equations, using 3 Si — A D 1 , AN — EJV coupled channels 

Chiral Quark Model local interactions, and ttN and coupled tvA — 7rE separable interactions fitted 
O ■ 

. to the position and decay parameters of the A(1232) and E(1385) resonances, respectively. The 

results exhibit a strong sensitivity to the p-wave pion-hyperon interaction, with a 7rAiV quasibound 
state persisting over a wide range of acceptable parametrizations. 
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I. INTRODUCTION 



The success of the nonrelativistic quark model (QM) during the 1970s in reproducing the SU(3) octet and decuplet 
baryon masses in terms of 3q configurations was followed by QM studies of 6q configurations that aimed particularly at 
elucidating the baryon-baryon short-range dynamics and making related predictions for dibaryon bound or quasibound 
states. It is remarkable that decades of experimental searches for dibaryons have so far yielded no unambiguous 
evidence for a dibaryon state. In the nonstrange sector, where the quark cluster calculations for L — 6q configurations 



m ■ n 

[l| suggest only a weakly bound AA dibaryon with (J P ,I) = (3 + ,0), there is a recent indication from np —> dmr 



reactions at CELSIUS- WASA for a resonance structure at Mr « 2.36 GeV and Tr m 80 MeV that might suggest a AA 
dibaryon bound by about 100 MeV, but still about 200 MeV above the dmr threshold [2J. In the strange sector, Jaffe's 
dibaryon H Q with strangeness S ~ —2 and quantum numbers (J P ,I) = (0 + ,0) which was predicted as a genuinely 
bound state well below the AA threshold, perhaps the most cited ever prediction made for any dibaryon, has not been 
confirmed experimentally to date in spite of several extensive searches [4] . Another equally ambitious early prediction 
was made by Goldman et al. [5|, also using a variant of the MIT bag model, for ( J p ; /) = (1 + , 2 + ; |) S = —3 dibaryons 
dominated by flN structure and lying below the SA threshold. More realistic quark cluster calculations by Oka et 
al. Q], applying resonating group techniques, did not confirm Jaffe's deeply bound H, placing it just below the 5 A 
threshold as a resonance about 26 MeV above the AA threshold. The underlying binding mechanism common to all 
of these orbital angular momentum L = configurations is the dominance of the color-magnetic interaction for gluon 
exchange between quarks, a feature emphasized by Oka [7( who systematically studied L = dibaryon configurations 
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that may benefit from a short-range attraction. Foliowing earlier quark cluster calculations [1J, |6j , these calculations 
resulted in no strange dibaryon bound states, and for the f2A-dominated S = —3 bound-state configurations predicted 
in Ref. [5], in particular, only a (J P ,I) = (2+, |) quasibound state resulted. 

For strangeness S = — 1, which is the focus of the present work, no L = dibaryons have been suggested for 
the lowest energy I = \ AN - EN coupled channels, where the long-range pion exchange interaction is dominant, 
particularly for the 3 S\ — 3 D\ system through the tensor component. Although old K d — >• 7r Ap data [8| had 
suggested resonant Ap structures at the UN threshold and 10 MeV above it, a (J P ,I) = (1 + , i) £A quasibound 
state is not necessarily required in order to reproduce the general shape of the Ap spectrum as shown by multichannel 
Faddeev calculations [9|,[l0|. Several low- lying L = 1 AN resonances were predicted in singlet and triplet configurations 

but negative results, particularly for the singlet resonance, were reported in 
13 1 near the £A threshold. At higher energies, Oka's analysis [7| drew 



in a QM study by Mulders et al. 
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dedicated K~ initiated experiments 
attention to a (J P ,I) = (2 + , |) dibaryon predominantly of a £(1385)A — £A(1232) coupled channels structure 
resonating about the £A(1232) threshold, approximately 100 MeV above the lower £(1385)A threshold. We note 
that these two channels are substantially higher in mass, by about 300 MeV, than the S = — 1 thresholds of AN and 
EJV. 

In a recent paper [14] we studied within three-body Faddeev calculations the possible existence of a 7rAA quasibound 
state, driven by the two-body ( 2S+1 Lj,I) = ( 2 Ps , §) nN resonance A(1232) and the ( 2 P|, 1) ttA resonance £(1385), 
for a AA ( 3 Si, |) configuration, all of which were represented by means of single-channel separable potentials. The 
coupling to the pionless SAT channel, with threshold about 60 MeV below the 7rAA threshold, was disregarded. It was 
felt that this coupling was mostly responsible for the width of the 7rAA quasibound state. The three-body channel 
(J P ,I) = (2 + ,|) was selected since all the angular momenta, spins, and isospins in this channel have maximum 
values and, therefore, it is likely to benefit from maximal attraction of both A(1232) and £(1385) resonances. This 
opportunity is unique to strange and charmed systems: a similar choice of ( J p , I) = (2 + , 2) for irNN, with each ttN 
pair interacting in the (| + , |) resonating channel, implies a (1 + , 1) Pauli-forbidden NN configuration. In terms of 
dibaryons, the ttAA (J p ,I) = (2+, §) quasibound state is a deeply bound £(1385)A - AA(1232) L = dibaryon, 
at energy considerably below the (J P ,I) = (2+, |) £(1385) N - £A(1232) L = dibary on suggested by the quark 
cluster model of Ref. [7| . 

Whereas the interactions in the pion-baryon resonating channels in first approximation are adequately represented 
by rank-one attractive separable potentials, the baryon-baryon interaction requires a rank-two separable potential to 
simulate both the attraction and repulsion that meson-exchange models normally yield. Indeed, we found a strong 
dependence of the calculated 7rAA binding energy on the balance between repulsion and attraction in the 3 <5i A A 
channel [ijj]. It is therefore suggestive to consider a more realistic hyperon-nucleon interaction in the J p = 1 + 
channel. In the present work we used the hyperon-nucleon (YN) Chiral Quark Model (CQM) interaction described in 
Refs. QQ in terms of 3 Si — 3 D\ , AA — £A coupled channels local potentials. For consistency, we also generalized 
our previous single-channel model of £(1385) as a nA resonance to a family of pion-hyperon (irY) interaction models, 
in terms of 7rA — 7r£ coupled channels separable interactions fitted to the position, width and decay branching ratios 
of £(1385). Furthermore, we studied the dependence of the calculated 7rAA binding energy on the irY interaction. 



In our previous work 14j, based on separable potentials, we considered both a nonrelativistic and a relativistic 



three-body formalism from which we deduced that the nonrelativistic results do not change much when the relativistic 
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formalism is used instead. This is relevant for the validity of the results of the present work which are based on the 
hyperon-nucleon interaction derived from the CQM within a nonrelativistic formalism. Therefore, in the present 
calculation we consider only a nonrelativistic framework. The results of the present three-body Faddeev calculations 
leave wide room for the existence of a (2+, |) wAN quasibound state indicating, however, a strong dependence on the 
short-range behavior of the least known irY and YN two-body subsystems. 

The plan of the paper is as follows. In Sec. II we describe the two-body interactions in the pion-nucleon, pion- 
hyperon and hyperon-nucleon subsystems. In Sec. Ill we derive the Faddeev equations of the pion-nucleon-hyperon 
system. Finally, we discuss our results in Sec. IV and summarize the work in Sec. V. 



II. THE TWO-BODY SUBSYSTEMS 



We will denote the hyperon, nucleon, and pion as particles 1, 2, and 3, respectively, and refer to the two-body 
subsystems by a subscript for the spectator particle. Thus, pion-nucleon is subsystem 1, pion- hyperon is subsystem 
2, and hyperon-nucleon is subsystem 3. The conventional reduced masses are given by 

= WjTOfc ^ a = m t (m 3 +m fc ) 

1 nrij + rrik ' 1 rrii + rrij + rrik ' 

where a superscript a = A, E has been added to indicate whether particle 1 is a A or a E hyperon and, obviously, 

Vi=Vi=Vi = : ■ (2) 

to?,- + rriN 



However, an average hyperon mass 



was used in the following reduced masses: 



m A + m E 

my = (3) 



TOA/(m 7r + my) m, mjv+my 

i/ 2 = , v z = . (4) 

mjv + mjr + my + + my 

The irY and YN amplitudes are given by 2 x 2 matrices, to account for the coupling between nA and 7rE and 
between AN and T,N, respectively. The irN amplitude in the three-body system is also given by 2 x 2 matrix, since 
the energy dependence of the two-body subsystem depends on whether the spectator particle is a A or a E. 



A. The pion-nucleon subsystem 



Since the irN subsystem is dominated by the A(1232) resonance, a rank-one separable interaction is considered 
sufficient: 

<pi\Vi\p' 1 >= im{pi)gi(Pi) , (5) 

so that the corresponding two-body t-matrix is given by 



<pi|ti(£)| Pl >=ffi(pi)7i(£)ffi(pi) , 



(6) 
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where E = Pq/2t)i with po the correct relativistic ttN center of mass (cm.) momentum and 

rrHE) = l/7i - r P ld Pl - 9 t^\. . (7) 
Jo E -p(/2r]i+ie 



The form factor gi(pi) was obtained from a very good fit of the P33 phase shift 17j for < 7i a b < 250 MeV in the 
form 



.91 



(p 1 )=p 1 [e-rtfi+A 1 p 2 1 e-^ a l] , (8) 



with 71 = -0.03317 fm 4 , A 1 = 0.2 fm 2 , /3 X = 1.31 far 1 , and a x = 3.2112 fitn -1 . 

In the three-body calculation, when the ttN subsystem is embedded in the irYN system, the energy argument 
of the isobar propagator ti(E) depends on whether the spectator hyperon is a A or a E, so that the separable irN 
ampitude ([6]) takes the form 

where 

ra(gi) = n(E-S aS AE-ql/2^); a = A, E, (10) 
with gi the relative momentum between the hyperon and the ttN subsystem and 

AE = TO S - TO A . (11) 



B. The pion- hyperon subsystem 

The ttY subsystem is dominated by the £(1385) p-wave resonance which decays mainly into 7rA and 7rE with 
branching ratios of (87.0±1.5)% and (11.7±1.5)%, respectively [l8]. To account for the coupling ttA~ 7rE, we assume 
a coupled channels separable interaction: 

< P2\vf\p' 2 >= 72 .9 2 > 2 ).9 2 V 2 ); a,0 = A, E, (12) 

so that the corresponding two-body t-matrix is given by 

< p 2 \t 2 (E)\p' 2 >=g%(p2)T 2 (E)g$(p' 2 ); a,/3 = A,E, (13) 

with 

7 E -p 2 , /2r)£ + te J E - AE - p 2 t /2r)% + it 

Again, _E = p\j2r\ 2 where po is the correct relativistic it A cm. momentum and A.E is chosen such that the it A 
momentum at the 7rE threshold has its correct value, that is 

AE = [(m s + to.) 2 - (toa + to.) 2 ][(to s + to.) 2 - (m A - to.) 2 ] 

877^(7715] + TO.) 2 
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TABLE I: Five choices (A-E) of form factor parameters for the coupled channels nY subsystem, Eqs. (|12[) and (|18[) . The last 
line lists a single-channel 7rA sixth model (F) with c-z = 0. The last column lists values of the r.m.s. momentum [see Eq. (|19|) ], 



Model 


A 2 (fm 2 ) 


«2 (fm 1 ) 


72 (fm 4 ) 


C2 


V<i>2 >S2 (fm" 1 ) 


A 


0.8 


2.41372 


-0 00604931 


890227 


4.11 


B 


1.0 


2.29039 


-0.00552272 


0.925591 


3.91 


C 


1.2 


2.20024 


-0.00501334 


0.956818 


3.76 


D 


1.5 


2.10192 


-0.00433208 


0.997300 


3.60 


E 


1.8 


2.03076 


-0.00375829 


1.03166 


3.48 


F 


3.21 


2.20024 


-0.00149204 





3.79 



The nY i-matrix (|T3j) in the ttYN system may be written in compact notation as a 2 x 2 matrix 

rjv(<fc)(<s£| < ff2 E |), (16) 



l.9 2 A > 



where 



T N {q 2 ) =T2{E-qt/2v 2 ). (17) 

The form factors g\ (p 2 ) of the separable irY p-w&ve potentials were taken in the form 

ff 2 A (p 2 ) = pa(l + A 2 p 2 2 )e-^ a l gf(p 2 ) = c 2 g^(p 2 ). (18) 

Solutions exist for all values of A 2 between and 00. Therefore, in order to fit the position, width and decay branching 
ratios of £(1385), we have at our disposal four free parameters: A 2 , a 2 , 72 and c 2 , which provide for varying one 
of these while adjusting the other three to the three pieces of data. We thus constructed five models (models A-E) 
by considering five values of the parameter A 2 , as shown in Table HI We also constructed a sixth model (model F) 
which shares the same range parameter a 2 with model C but which neglects the coupling to the 7rE channel {c 2 = 0), 
as was done in our previous calculation [14]. It is instructive to classify the various ttY interaction form factors 
g\ ' (p 2 ) according to their root-mean-square (r.m.s.) momentum, using the following expression for the mean-square 
momentum < p\ > ff2 : 



2 _ Jo 92(P2) v\ d 3 P2 _ 2 Maj , g 2 
Jo .92(^2) d 6 p 2 A 2 a^ + 2 

where the approximation owes to 2A 2 a 2 >> 1. The resulting values of the r.m.s. momentum, listed in the last column 
of TableCD are close to v / < v\ > 9 2 ~ 3-8 fm" 1 « 750 MeV/c. For comparison, ^< p{ > gi = 5.55 far 1 w 1100 MeV/c 
for the irN form factor gi(pi) of Eq. ([8]). 1 



1 This high-momentum value for g\ does not rule out a spatial size of order 1 fm for A(1232). Indeed, if gi(r\) is the Fourier transform 
of gi (pi), for 1 = 1, then J< r\ > gl = 0.875 fm. 
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C. The hyperon-nucleon subsystem 

The YN interaction derived from the chiral quark model is a local potential obtained by application of the Born- 
Oppenheimer approximation to the chiral quark-quark interaction (consisting of confinement, one-gluon exchange, 
5seudovector-meson exchange, and scalar-meson exchange) with a fully antisymmetrized six-quark wave function 
In the case of the J p — 1 + , I = \ channel, it leads to the following system of coupled equations: 



15 



pOO 

tu"(&i,P3\ E )= Vu"(P3,P3)+ EE/ Pz d Pz V M>(P3,Pz) 

■v— A y. Qi— n 9^0 

-t],%(p 3 ,p s ';E): a,/3 = A,E, (20) 



7=A,S^'=0,2 ' 
1 



E - S^AE - p' 3 '/2r)2 + ie 

with a, ft = A, E, 1,1" = 0,2 and E = P q/2?] 3 , where po is the correct relativistic AiV cm. momentum, and AE is 
chosen such that the AN momentum at the E7V threshold has its correct value, that is 



x r [(mg + m N ) 2 - (m A + m JV ) 2 ][(m s + m N ) 2 - (m A - m N f] 

8773 (ms + m N y 



The YN i-matrix (|2T)|) may be written in compact notation as a 2 x 2 matrix 



^AA ^AE 



where each YN i-matrix t a ° includes, in addition, a coupling between S (£ = 0) and D (t = 2) waves. 

III. THE THREE-BODY EQUATIONS 

The Faddeev equations for the bound-state problem 

T^J^tiGoTf, 1,7 = 1,2,3, (23) 

couple the amplitudes Ti, T 2 and T3 together. Eliminating the amplitude T 3 in favor of T\ and Ta> one obtains 

Ti = iiGoiaGoTi + (ii + iiG i 3 )G T 2 , (24) 

T 2 = t 2 G i 3 G T 2 + (t 2 + t 2 Got3)G Ti, (25) 
where, in order to allow for the Y = (A, E) specification, one has 

( G o \ 

Go = I , (26) 
V G o) 

Since the two-body amplitudes t\ and ti are separable [see Eqs. (HJ) and tfl6])] . the three-body amplitudes T\ and T 2 
are of the form 

T i = l.9i > I XA I r>7) 



X 



T 2 \\ 9 i > \X,. (2.S) 

where the subscript of the amplitude X indicates which particle is the spectator. Substitution of ([27]) and ([28]) into 
(l24l) and (ESI) leads to 

X A \ /r A o\ Vg a ; aa G a G^gA > IX A 
XJ rj 511 \G ^ SA G A G$t™G$J \X 1 







: 









+ G A i AA G A G A t AS G s W|<7 A > 
Gq i SA G A G S + G^G S ] Lf > 



AV, (29) 



v l A| SI 1 I + 1 I . / A A 

A W = 7» ( < ff A| < gf\ ) ( Mfll > 



/ x, G A t AA G A G a £ as Gq \ / | 5 A > i 

( < 9 2 A | < 5 2 S |) 30 
* J ' G$t Gq G^ G z \\ g f > 



which are integral equations in one continuous variable given explicitly by 

p oo 

x a (qi)= r a ( qi ) J2 / glWiff^fa.gi^Wi) 

/•OO 

+r a (qi) q 2 2 dq 2 K aN ( qi ,q 2 )X N (q 2 ); a = A, E, 
Jo 



(31) 



/"OO 

^Arfe)= rjvfe) ^ / q 2 1 dq 1 K Na (q 2> q 1 )X a (q 1 ) 

a=A,E- 7 ° 

/•OO 

+Tiv(S2) / q 2 2 dq 2 K NN (q 2 ,q' 2 )X N (q' 2 ). (32) 
Jo 

The kernels of these integral equations are given in the Appendix. 



IV. RESULTS AND DISCUSSION 



We applied the formalism of the previous section, using six different versions of the YN interaction obtained from 
the CQM, all of which reproduce equally well the experimental low-energy YN data [3, Q|. Results are listed 
in Table HI] from where it is clear that the ttAN binding energies are substantial for nY models with A 2 < 1 fm 2 . 
Generally, the higher the r.m.s. momentum of the ttY form factor g 2 , the stronger is the binding, as demonstrated 
in the table. Irrespective of which ttY model is chosen, the YN interaction always produces repulsion, thus lowering 
the calculated binding energy, as demonstrated by the results listed in the last line which corresponds to switching 
off the YN interaction. This repulsive YN effect owes its origin to the high-momentum components of the ttB form 
factors which within the three-body calculation highlight the short-range repulsive region of the YN interaction. 

To demonstrate the model dependence of the three-body calculation within a given ttY model, we assembled in 
Table [Til] several binding energy results based on model C and also on its limitation to the ttA channel (model F 
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TABLE II: Binding energy of 7rAiV (in MeV) for five irY interaction models (A-E of Table HJ and six CQM versions of the 
3 Si — [i Di YN interaction fitted to given AiV scattering length a and effective range ro (both in fm). The momentum pi a b(<5 = 0) 
is the A laboratory momentum (in MeV/c) where the 3 Si AN phase shift changes sign. The last line corresponds to switching 
off the YN interaction. 



a 


ro 


Plab(<5 = 0) 


Model A 


Model B 


Model C 


Model D 


Model E 


-1.35 


3.39 


987 


147 


99 


65 


30 


6 


-1.40 


3.32 


1011 


147 


99 


66 


30 


6 


-1.64 


3.09 


1146 


150 


102 


68 


32 


8 


-1.71 


3.03 


1198 


150 


102 


68 


33 


9 


-1.78 


2.98 


1272 


151 


103 


69 


33 


9 


-1.86 


2.93 


1446 


152 


104 


69 


34 


10 








170 


120 


84 


47 


21 



TABLE III: Comparison of irAN binding energies (in MeV) calculated within the nY model C and its irA limit model F (see 
Table |T| for YN coupled channels and AN single-channel models with a — —1.40 fm, and for no YN interaction. 



ttA, AiV 


ttA, YN 


irY, no YN 


tyY, AN 


ttY, YN 


93 


96 


84 


73 


66 



listed in Table [I]). The YN models included in this table invariably give a — —1.40 fm, whether limited to the AiV 
3 Si single channel or extended to the 3 Si — 3 Z?i AiV — EiV coupled channels. Comparing the first two entries to 
each other, we conclude that the extension from a single AiV channel to AiV — EiV coupled channels has very little 
effect (about 3 MeV additional attraction) within the ttA model F. In contrast, for the full irY model C, as in the last 
two entries, the extension from AiV to YN models has a somewhat larger effect (about 7 MeV repulsion) and in the 
opposite direction. Within the ttY model C, the full coupled channels YN interaction contributes 18 MeV repulsion 
(third and fifth entries in Table |TTT| to the three-body binding energy. Similar results hold for all other irY models. 

To discuss the model dependence of the three-body calculation within a given YN CQM, we follow Ref. in 
singling out pinb(S — 0), the momentum where the AiV 3 Si phase shift changes sign from attraction outside to 
repulsion inside, as a measure of the repulsive YN effect. We notice in Table [TT1 that the CQM values of pi a t>(<5 = 0) 
are considerably larger than those obtained by other models 20j-|22|, signifying less repulsion in the CQM. In order to 



test whether the apparent lack of repulsion in the CQM YN interaction is responsible for the large binding energies 
obtained for A2 < 1 fm 2 , we added to the CQM with a = —1.40 fm and ro = 3.32 fm a short-range potential in the 
3 Si AiV partial wave of the form 

V(r) = lR 1A , (33) 

r r 

with /3r — 10 fm -1 and 7^ > 1000 MeV fm, while the attractive term was adjusted to maintain the AiV scattering 
length a = —1.35 fm and the effective range ro as close as possible to 3.39 fm, so that the YN observables are not 
changed noticeably. The overall effect of V(r) is repulsive, as demonstrated in Table ITVl for the nY model A, with 
A2 = 0.8 fm 2 , where it is clearly seen that increase in the strength of the repulsive term lowers the value of pi a b(^ = 0) 



as well as lowering the 7rAiV binding energy. 
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TABLE IV: Binding energy of tyAN (in MeV) for the ttY model A (A 2 = 0.8 fm 2 ) and the CQM YN interaction plus a 
short-range AN potential V(r), Eq. (|33[) . with scattering length a — —1.40 fm and effective range ro = 3.32 fm. The strength 
parameters 7 are in units of MeV fm, the inverse range parameter /3a is in units of fm -1 , /3r = 10 fm -1 , and pi a b(<5 = 0) is the 



laboratory momentum (in MeV/c) where the 3 Si AN phase shift changes sign. 

JR JA Pa Pi^{8 = 0) B{A 2 = 0.8) 

1000 240 5.371 873 107 

2000 530 5.811 846 88 

3000 720 5.749 822 70 

4000 990 5.928 810 59 

5000 1260 6.056 802 51 

6000 1360 5.921 788 39 

7000 1670 6.086 775 34 



V. SUMMARY 

In this work, we have extended the Faddeev equations study of a (J p ,1) — (2 + , |) quasibound ttAA^ state |14| from 
a 3 Si AN single-channel to 3 Si — 3 Di , AN — UN coupled channels, and from a irA single-channel description of 
£(1385) to 7rA — 7rE coupled channels description. Local interaction potentials given by the CQM were used in the YN 
sector, whereas one-rank separable potentials were used in the ttB sectors. We have shown within a nonrelativistic 
version of the Faddeev equations, but using semirelativistic kinematics, that the ttAN system is bound under a wide 
choice of parametrizations of the irY interaction form factor. The form factors of the ttB subsystems are sufficiently 
short ranged such that the pion undergoes almost coherently attraction to both baryons. The short-ranged repulsion 
between the two baryons in the CQM is insufficient to overcome the attraction gained by the pion unless the CQM 
is modified arbitrarily at very short distances to do this job. Altogether, the acceptable model dependence of the irY 
interaction form factor, and the uncertainty of the short-range behavior of the YN interaction, leave plenty of room, 
theoretically, for a quasibound S — — 1, (J P ,I) = (2+, |), 7rA7V dibaryon. 

Before closing we list several production reactions, where the first two were already discussed in our previous paper 
[lJ ]. in which to search for this S = — 1 dibaryon here denoted V: 

K-+d^V~+n + , n-+d^V~+K + , (34) 

p + p -> V + + K + . (35) 

Correlated with the missing mass spectrum of the V dibaryon, for a forward outgoing meson, one should look for 
YiN decays that can be assigned to a UN resonance with invariant mass M-p- Total cross sections for the associated 
strangeness production pp — > T,NK + near the hyperon production threshold have been reported from Juelich, for E°p 



by the COSY-11 Collaboration [23J, for E+n, also by COSY-11 [2J, and by the ANKE Collaboration |25j and the 
HIRES Collaboration [26 1 , with conflicting results among all these E + n reports. Old DISTO data for the reaction 
pp — > ApK + have been analysed to search for an intermediate K~pp quasibound state, with the astounding report of 
a broad resonance at the nYN threshold 27]. Of course, this I = \ resonance cannot be assigned to a I = | ttAA'' 
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quasibound state, but forthcoming data from the FOPI detector Collaboration at GSI [28] could be analysed also 
with respect to a SiV rather than a Ap final state. 
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Appendix: Expressions for the kernels of the integral equations Eqs. (|3ip and (|32p 

We provide here detailed expressions for the kernels appearing in the integral equations Eqs. (|31[) and 

l r i 



K^{q u q' x ) = I Y, r<& d( &f dcos0 I dcose ' 

D 01 — n n " J 1 J 1 



W=0,2 ' 

lb AlOil 



Xffi(pf )G?(pf , qimA^Uqu q 3 , co S 6) 
xt^^E-qi^b^ (36) 



K aN ( qi ,q 2 ) = \ J dco S eg 1 (p°)G°(p°,q 1 )b? 2 A{° 2 Z(quq2,cos6)g%(p« 
+ J ^ ^2 J <£ d< & J i dense dcosO' 



x g x (pf)G% (p? , qi)b? 3 A™% ( Ql , q 3 , cosfl) 

X$?(p£,j/f ; E - g 3 2 /2^ 3 )44' 2 n i° 92, cos0')G£(p£, q 2 )gl(pl), (37) 



K Na (q 2l q 1 ) = K aN (q 1 ,q 2 ), (38) 



pOO pi pi 

K NN (q 2 ,q' 2 ) = ~ E E / ^ I dc ° s9 I dc0s9 ' 



' £,£'=0,2 a,/3=A,E 



J-l J-l 



xtf (P?)Go (Pa, q2)b^AH%{q 2 , q 3 , cosfl) 

xt^(p«,ptE - ql/2v 3 )bl 2 AQ™( q3 , q' 2 , cos9')G^p' P 2 , q' 2 )g$(p' 2 ), (39) 



with 



G ^= E-S aS AE-p 1 /2r 1? ^/2^ + ^ * = A < S " ™ 

£ X £ X £ ■ X • £ X 

The orbital angular momentum recoupling coefficients AA £ Q J J (qi , qj , cos(9) = A.? ' (qj,qi,cos9), and isospin 
recoupling coefficients = are calculated by consideration of a cyclic pair ij. (The spin recoupling coefficients 
are all equal to 1.) For isospin we have 
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b% = (2z, + l)(2ij + 1) ^(r^/ri; Uij), (41) 

where FT is a Racah coefficient. If a = A then ti = 0, T2 = ^, T3 = 1, i\ = |, 12 = 1, 43 = 5, and / = | so that 
= 6^ = = 1. If a = S we have instead n = 1 so that 6f 2 = \/5/6, 6f x = — v/5/3, and &2 3 = — l/y/6. 
The orbital angular momentum recoupling coefficients are given by 

xT e . mj T x . M _ m .cos(-M6-m i 9f + m j 9f); a = A,E, (42) 
where Te m — for odd values of £ — m, and 



_ (-l)('+ m >/y (21 + l)(l + m)\(£ - m)\ 
lm ~ 2 e [(l + m)/2}\[(e-m)/2}\ : 

for even values of I — m. The angles Of and 8" are obtained from 



qjCOsO + qta^ 



(43) 



cosC = ^ -^-!L, (44) 



Qi COS0 + QjO.% , s 

cos6 f = ± ^Ul^ 45 

Pi 



P? = ^q 2 3 +{q l af J Y+2q t q 3 af J coBe, (46) 
Pj = ^<lf + (<lja«y+2q l q J a«cose, (47) 

where 



(48) 



ti 1 ' it 1 1 

J rrij + ?7ife J mi + TOfe 

with mi = m Q ; a = A, S. Equations ()46[) and (|47jl provide also the relative momenta appearing in Eqs. (|36p ~([39 
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